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We apply negative dimensional integration method (NDIM) to three out- 
standing gauges: Feynman, light-cone and Coulomb gauges. Our aim is to 
show that NDIM is a very suitable technique to deal with loop integrals, 
being them originated from any gauge choice. In Feynman gauge we per- 
form scalar two-loop four-point massless integrals; in the light-cone gauge 
we calculate scalar two-loop integrals contributing for two-point functions 
without any kind of prescriptions, since NDIM can abandon such devices - 
this calculation is the first test of our prescriptionless method beyond one- 
loop order; finally, for the Coulomb gauge we consider a four propagator 
massless loop integral, in the split dimensional regularization context. 
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1. INTRODUCTION 

Perturbative approach for Quantum Field Theory in any gauge deals with Feyn- 
man diagrams, which are expressed as Z)-dimensional integrals. The success of such 
approach can be understood from the comparison between the a = ^(17 — 2) measure 
for the electron, 

QThe = 1159652201. 2(2. 1)(27.1) x lO^^^ 

asxp = 1159652188.4(4.3) x 10"l^ (1) 



see for instance [1]. 
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This is the best motivation for studying Quantum Field Theory, no physical 
theory can give such accuracy in any measurement. In other words, it is the very 
best we have. 

In section 2 we discuss some 2-loop 4-point functions, namely, on-shell double 
boxes with 5 and 6 massless propagators; section 3 is devoted to non-covariant 
gauges: the light-cone and Coulomb ones. The integrals we study for the former 
have 7 propagators (2-loops) and the latter is 1-loop and have 4, however it is also 
complicated since we have to use split dimensional regularization(SDR). In the final 
section, 4, we present our concluding remarks. 

2. FEYNMAN GAUGE: SCALAR TWO-LOOP FOUR-POINT 
MASSLESS INTEGRALS 

Of course, covariant gauges are the most popular, in what we could call "gauge 
market" [2] . Several methods were and are still developed to evaluate complicated 
Feynman loop integrals, being them concerned with analytic or numerical results[l, 
3, 4], all in the context of dimensional rcgularization[5]. 

Our work in concerned with the application of negative-dimensional integration 
method (NDIM). It is a technique which can be applied to any gauge, covariant or 
non-covariant alike. The results are always expressed as hypcrgeomctric scries which 
have definite regions of convergence allowing one to study the referred diagrams or 
process in specific kinematical regions of external momenta and/or masses. 

On the other side, NDIM has a drawback: the amazing number of series - in the 
case where one is considering massless diagrams - which miist be summed. When 
such sums are of gaussian type, it is quite easy to write a small computer program 
that can do the job algebraically. However, when the series are of superior order, 
p+iFp, for p > 2, there are no known formulas which can reduce it to a product of 
gamma functions for any value of its parameters. Despite this technical problem, 
NDIM proved to be an excellent method[6, 7, 8, 9]. 

A question which is often arised is: what is more difficult to handle, graphs with 
more loops or graphs with more legs? In our point of view, i.e., in the context of 
NDIM, the greater the number of loops the heftier the calculations will be needed 
to solve it. We will consider in this section a diagram which has both (great number 
of legs and loops, four and two respectively), a scalar two-loop double-box integral 
where all the particles are massless and the external legs are on-shell. 

2.1. Double box with 5 and 6 propagators 

Let us consider the diagram of figure 1. Consider as the generating functional 
for our negative- dimensional integral the gaussian one, where all external legs are 
on-shell. 



Gb = d^q d^r exp [-aq^ - (3{q - pf - j{q - p - p'f -e{q-r- pif - (tyr"^ 



■u{q-rf] 



(2) 




(3) 
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FIG. 1. Double-box with six propagators. 



where (s, t) are the usual Mandelstam variables and we use s + t + u — 0. Observe 
that in the particular case where a = we recover the gaussian integral for the 
diagram of figure 2. We also define A = a9 + acj) + aiv + PO + (3(j) + f3uj + j6 + -f(f) + 
juj + (jxu + 94>. 

The usual technique reveals that there are thirteen sums and seven equations. 
From the combinatorics one can solve such constraints in f7f6 different ways. Of 
course several systems have no solution - not even in the homogeneous case - and 
from our previous works we know that some results are rt-fold degenerated and 
others are related by analytic continuation. The result for the integral in question 



BOX= / d^qd^r {q^y{q-pf^{q-p-p'f^{q-r-pif\r'')"'{q-rf'^,{A) 

Xi\X2\Yi\...Yg\Zi\Z2\ 



{-T:)°i\mi\m\n\T(l - - D/2) V ■ — ■ — ■ — -,5, (5) 

aU=0 " 



where "all" means {Xi, X2, Fi, Fa, 1^3, n, F5, Y^.Yj, Y^, Yq, Zi, Z2} and 

BOX = BOX{iJ,k,l,m,n), 

must be understood and delta represents the system of constraints. The above 
expression can be expressed, in principle, as a seven-fold hypergeometric series, 
there are three possibilities, 

Ti...\z,z-\1), .^(...^,1), and T{...\z-\1), (6) 



where z = ~s/t. Some series with unit argument, if they were gaussian can be 
summed up. However, a hypergeometric function is meaningful only if the series 
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which defines it was convergent. Since the first possibihty cannot be convergent we 
disregard it. 

Among the 624 total sohitions of the system of constraints we look for the simplest 
solution, namely, the one in which we can sum the great number of series. It is not 
diflacult to find it using computer facilities, 

BOX^^{i,j,k,l,m,n) = h{i,j,k,l,m,n) sF2{{l}\z), (7) 

where the five parameters are quoted in the table, at = i+ j + k + l + m + n + D 
and 

h{i,j,k,l,m,n) = TT''{t^r'{-j\ab){-l\ab){<Jb + D/2\-2ai,~D/2) 

X {-m\l + m + n + D/2){i + j + k + m + D\ -m-D/2) 
X (Z + m + n + £i| -/-n-D/2), (8) 

besides this one, wc can have Appel's, Lauricella's and even more complicated 
hypergeometric functions. Moreover, 

{x\y) = {x)y = ^^^^ . 

If we remember that the final result should be the sum of linearly independent 
serics[7, 10], we can rightfully ask if one is not missing two other 3_F2 functions. 
According to Luke and Slater [12], the difi'erential equation for pFq has p linearly 
independent solutions, so we should write a sum of three terms. On the other hand, 
according to Njzirlund[13], if the difference between an upper parameter and a lower 
one was an integer number, then some series do no exist — we used this theorem in 
[7]. So, eq.(7) is the final result for the referred integral in the region where \z\ < 1. 
The expression for the same graph outside this region can be obtained making the 
substitutions, 

s <-> j ^ k, I n, (9) 

so we have other 3F2 hypergeometric function as the result for |2;| > 1. 

TABLE 1 



Pcirameters of hypergeometric functions 3F2 representing box integrals 



Parameters 


3F2({1}|^) 


zF^{{2}\z) 


a 


-k 


-k 


b 


—n 


—n 


c 


-0-6 


t 


e 


1 + j - 




f 







LOOP INTEGRALS IN THREE OUTSTANDING GAUGES 



5 



Another solution for the Feynman integral can be written as a triple hypergeo- 
metric series, 



{i + k + l + m + D\Y45 + 2Ye){i + l + m + n + D\Y45e)z^^^<' 



X 



(1 - j + l\Y56){i +j + k + m + D\Y46){1 + m + n + £>|F456) 



where 



/g = (-l\j)(l + m + n + D\i){ab + D/2\-j-n-D/2){-n\l + 2n + D/2) 

x{i + j + k + m + D\ - m - D/2){-k\j + k- ab){-m\2m + D/2), (11) 

obviously the above series converges if \z\ < 1, besides other possible condition 
on z. So there is an overlapping between the regions of convergence of BOX and 
BOX 3, so there exists an analytic continuation formula which relates both. As fas 
as we know textbooks do not show formulas relating triple hypergeometric series 
with simple ones. 

We have also 4- fold series, 



BOXi = 



^DA^a,-if V {-m2A7){m + D/2\Y,,j){-i\Y,) 

■'^^^ ri!y4!W2!(l +^6 - Z|yi247)(l + J - 1\Yi27) ^ ' 

{i + l + m + n + D\Y4r){i +j + m + n + -P| Y124 + 2Yr)z^^^''-' 
{i + l + m + n + D\Y4r)il + m + n + D\Yi4r) 



where 



/4 = {-j\l)il + m, + n + D\i + j-l){ab + D/2\-k-l-D/2){-m\2m + D/2) 
x{-k\k + l-ab){-n\l + 2n + D/2). (13) 

Observe that the two previous results are singular when j — I = integer, since 
we have T{j — I) or T{1 — j) in the numerator. However, such singularity cancels 
if one consider propagators exponents in the analytic regularization context, i.e., 
introducc[7, 11] for instance j = —1 + S, then expand the whole expression around 
6 = 0. Proceeding in this way the pole in 6 cancels. 

We have above reduction formulas which transform a hypergeometric function 
defined by triple and 4-fold series in a simpler function defined by a unique sum. 
These formulas are not in the textbooks on the subject. It is an original result. 

2.1.1. Double box with 5 propagators 

The graph of figure 2, is a special case of the previous one. In the gaussian 
integral a must be zero, so in the final result we must merely take i = 0, 



BOX^'' {0,j, k, I, m, n) = /(O, j, k, I, m, n) 3i^2({2}|z), (14) 
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FIG. 2. Double-box with five propagators. 

where the parameters are hsted in the table and we define cr^ — j + k + l + m + n + D 
and 

f{0,j,k,l,m,n) = n''{tY^{-j\<7l){-l\a',){<7l + D/2\-2a',-D/2) 

x{-m\l + m + n + D/2){j + k + m + D\-m- D/2) 
x{l + m + n + D\-l-n~ D/2), (15) 

we can proceed with the same substitutions (9) to obtain the result outside the 
region \z\ < 1. 

Finally, when all the exponents are equal to minus one, the 3F2 colapses to a 2F1 
which can be written as an elementary function. 

The results BOX 3 and BOX 4, in the same special case {i = 0), are hypergeo- 
metric series representations for the integral in question. Since they are different 
and depend on the same variable we must sum them in order to get a triple series 
representation for BOX^^O, j, k, I, m, n). 

3. NON-COVARIANT GAUGES: LIGHT-CONE AND COULOMB 

Recently there have been many works on non-covariant gauges, namely, light- 
cone[14, 15], Coulomb[16] and radial and axial gauges[17]. Despite they are not so 
popular as covariant ones, they have some important features which can help our 
study on certain physical problems. 

Light-cone gauge, as far as we know, is the only one where certain supersymmetric 
theories can be shown to be UV finite and possess a local Nicolai map [18]. Moreover, 
ghosts decouple from physical particles and one is left with a reduced number 
of diagrams. On the other hand, the price to pay seemed to be so high, since 
the gauge boson propagator did generate spurious poles in physical amplitudes. 
This problem was overcame when andelstam and Leibbrandt[19] introduced causal 
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prescriptions to treat such poles (there are also other causal prescription which can 
be implemented, proposed by Pimentel and Suzuki[20], known as causal Cauchy 

principal value prescription.) However, the famous ML-prescription necessarily 
forces one to use partial fractioning tricks and integration over components, which 
turn the calculations rather involved [21]. 

Negative-dimensional approach can avoid at all the use of prescriptions and pro- 
vide physically acceptable results, i.e., causality preserving results. The calculation 
we will present is the very first test beyond 1-loop order without invoking ML- 
prescription, as we called in [8] NDIM is a prescriptionless method. Still, integration 
over componenets and partial fractioning tricks can be completely abandoned as 
well as parametric integrals. The important point to note [8] is that the dual light- 
like 4- vector n* is necessary in order to span the needed 4-dimensional space[18, 22]. 
This is the very reason why our calculation for one-degree of covariance violation 
failed [6]. 

The second non-covariant gauge we deal with in this paper is the Coulomb gauge. 
Potential between quarks and studies on confinement are easily performed in this 
gauge[15, 16]. Besides, the ghost propagator has no pole in this gauge! As light- 
cone gauge. Coulomb also have problems with gauge boson propagator. In the 
former, loop integrals generated aditional poles; in the latter, such integrals are not 
even defined [23] since they have the form, 

J (16) 

such objects are the so-called energy-integrals. Doust and Taylor [23] presented a 
solution for this issue in a form of a interpolating gauge (between Feynman and 
Coulomb). Leibbrandt and co-workers [24] presented also a solution, a procedure 
they called split dimensional regularization, (SDR), which introduces two regulating 
parameters, one for the energy component and another for the 3-momentum one. 
So, the measure becomes, 

A = dq^d'^-^ "-^ d^q = dPq^d'^q, (17) 

in Euclidean space. 

NDIM can also deal with Coulomb gauge loop integrals, but it needs to make 
use of SDR. In this work we propose to apply NDIM to scalar integrals with four 

massless propagators. Our results are given in terms of hypergeomctric series in- 
volving external momenta, exponents of propagators and regulating parameters u> 
and p. 

3.1. The Light-Cone Gauge 

So far, we have tested our NDIM for integrals pertaining to one-loop class. Now 
we apply such technology to some massless two-loop integrals. Let us consider an 
integral studied by Leibbrandt and Nyeo[21], since they did not present the full 
result for it. 
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where in their calculation ML-prescription must be understood. On the other hand, 
in the NDIM context the key point is to introduce the dual vector n* in order to span 

the needed spacc[8, 18, 22]. If wc do not consider it, our result will violate causality, 
giving the Cauchy principal value of the integral in question, as we conclude in [6] . 
NDIM can consider lots of integrals in a single calculation. Our aim to is perform, 



Ar= / d''qd''ki{kly{qyiq-kif''iki-pf{ki-nr\q-ny{ki-n*Y, (19) 



we will carry out this integral and then present results for special cases, including 
Leibbrandt and Nyeo's C3, where i = —1, r = and the other exponents equal to 
minus one. Observe that the integral must be considered as function of external 
momentum, exponents of propagators and dimension, 

JV = J\f{i,j,k,l,m,r,s;P,D), (20) 

where P represents (p^,p"'",p~, ^(n • n*)), and we adopt the usual notation of light- 
cone gauge [2]. 

Our starting point is the generating function for our negative-dimensional inte- 
grals, 

Gn = j d^q d^k cxp [-afc2 - Pq^ - j{q - kf - 0{k - pf - (j){k ■ n) - uj{q ■ n) 
-ri{k ■ n*)] , (21) 

then after a little bit of algebra we integrate it, 

1 



Gn 
where 





D/2 






1 exp- 


.A 



-ffi/ - 92{p ■ n) - g3{p ■ n*) + g^i^n ■ n*) 



(22) 



92 

gi = {aj3 + a'y + j3-f)e, g2 = {(3(l> + ^u) + ^(j))9, g^ = {(} + -f)T]9, g4 = r]—, 

and X = aP + a^ + 0^ + pe + ^6. 

Taylor expanding the exponentials one obtain, 

N = {-.-W.kmm\rlsni-On-Dl2) ^ . . . X^Y^W^Y^l 

all=0 

z,\...z,\ 1, ir) ' ^'"^^ 

where (7„ = i + j + k + l + m + r + s + D and 5 represents the system of constraints 
(8 X 16) for the negative-dimensional integral. In the end of the day we have 12870 
possible solutions for such system! Most of them, 9142, have no solution while 
3728 present solutions which can be written as hypergeometric series. Of course 
several of these will provide the same series representation, these solutions we call 
degenerate. 
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We present a result for the referred integral as a double hypergeometric series, 

{an + D/2\Zi5){i + j + k + m + s + D\Z^^){D/2 + k\Zi) 



^ _^ Z^\Z^\{1 + i + j+k + an+ D\Z^5){j + k + s + D\Zi5) 

^, (.7 + s + D/2\Z5){i +j + k + r + DIZ45) f p^n ■ n* \ 

^ {l + i+j + k + m + r + s + D\Z45) \2p+p-) ' ^ ' 



where 



fn = {-m\ - s)i-i - j - k - D/2\ ^ an - D/2)ij + k + s + D\i ~ s + r) 
+ I + D/2){-k\ - J - D/2){-m\j +m + s + D/2) 
{l + r\— i — j — k — m — r — s — D) 
x(— j| — i — k — m — r — s — D), 

are the Pochhammer symbols and 

P„ = (p2)a„+i+j+fc+D(p+)i+m+.-a„(p-)!+r-a„ " \ (26) 



Now wc can consider the special case {i = l,j = k = I = m = s = —l,r = Q), 
studied in [21], 



^ n r(5 - 2D)T{D - l)r(D/2 - l)r(2 - D/2)T{D/2 - 2) 

r(l - D/2)T{D - 3) ' 

'n-n*\"-^ ~ {3D/2-4\Z45) 



\ / 2 * \ -^45 

r \ I 71 n • n \ 

(27) 



{D/2 - l\Zi){D/2 ~ 2\Z^){D - IIZ45) ( p'^n ■ n*^ 



{D - 2IZ45) V P+P- 

observe that it exibits a double pole, as stated by Leibbrandt and Nyeo[21]. 

3.2. The Coulomb Gauge 

We will present the full calculation of an integral which has four propagators, 

J(*,j, fc,m) = j d^q {q^r{q-pf'<l^\ci + pf"^, (28) 

in order to regulate the possible divergences originated by the energy component, 
SDR must be understood, namely, 

d'^q = dPqid^ci, (29) 

where D = p + w. 

The generating functional for our negative- dimensional integrals is the gaussian- 
like integral. 



Ge = y A exp [-aq^ - (3{q + pf - - 0(q + p)^] , 



(30) 
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which can be easily integrated, 

/ nR ^\ r (n -I- ^^(R -\- _1 

(31) 



(a + 7)(/3+^ 2' 



There are results given by double, triple, 4-fold and 5-fold hypergeometric series in 
the variable p^/p4 or its inverse. 

We will present two of such hypergeometric series representations, the first one 
is a 4-fold series, 

J.U 1km)- C.U lkm)Y^ H\X,.u){j+m + D/2\X,i) /p^N^^ 
Mt,J,k,m) - C4{i,j,k,m) 2^ Xi!X2!X3!X4! \pl) 

— 

^, + 7 + ni + ij/2\X,, - X,2'){-ni\X2) 

^ {l+j + k + m + D/2\X34:) 

{-j - p/2\Xi ~ X^){k + L^/2|Xi24) 

(1 - i - p/2|Xi24) 

j,k^ m), (32) 

where 

C4{i,j,k,m) = n''^\ply{p''yc-i^_j\_p/2){-j-m-p/2\-k-u;/2) 

X {-k\2k + u)/2){j + k + m + D/2 + u)/2\ -k- w/2), (33) 

where ac = i+ j + k + m + D/2 and the second a 5-fold hypergeometric series, 

T.a i k m) - CM i k m)Y^ j-i - j - p/2\2Xi + X2345) 
Mt,j,k,m) - C,{z,j,k,m) 2^ x,\X2\Xs\X4\X,\ 



pf/ {l + k + m + Lj/2\Xi45) 



(-l)^^(m + w/2|Xi345) 



ik + u/2\Xi245) 



(1 - j - p/2|Xi35)(l - i - p/2|Xi24) ' 

where 



(34) 



C,{i,j,k,m) = 7r^/2(p2)»+^-+p/2(p2)fe+™+-/2(_,|, + p/2){-j\i + j + p/2) 
x{-k\k + m + u /2){-m\k + m + w /2){i + j + p\ -a^- p/2) 
x{k + m + uj\-ac-uj/2), (35) 

observe that the above result is also symmetric in [i ^ j,k ^ m), which means in 
the loop integral, q^^ +p^. 

Another important point to observe is that the final result must a sum of linearly 
independent hypergeometric series[6, 7]. The above 5-fold series, J5, appears only 
one time whereas J4 is degenerate since several systems give its two hypergeometric 
functions. This must be considered if one wants to apply NDIM to more complicated 
diagrams which can in principle generate hypergeometric series representations even 
more involved. 



LOOP INTEGRALS IN THREE OUTSTANDING GAUGES 



11 



Moreover, the above expressions, J4 and J5, are related by direct analytic contin- 
uation, since both are convergent for Ip^/pH < 1. When one is considering simple 

hypergeometric function, several formulas are known; on the other hand, for rather 
complicated hypergeometric series, as we obtained — four and five-fold series — , 
there are very few of such formulas. NDIM can fill this gap, since it is the only 
method which provides hypergeometric series representations for Feynman loop in- 
tegrals, in different kinematical regions, and related by analytic continuation direct 
or indirect alike. 

The above hypergeometric series (only the series, not the factors!), J4 and J5, can 
be written as generalized hypergeometric functions [12] of four and five variables. 



T 



6:0;0;0;0 
2:0;0;0;0 



{-% ■ 1, 1, 1, 1), (i + m + D/2 : 0, 0, 1, 1)(1 + j + m + p/2 : -1, -1, 1, 0) 
(l+i + fc + m + £>/2 : 0,0,1,1) 

(-m : 0,1,0,0) 
(l-z-p/2: 1,1,0,1) 



(36) 



and 



{-i-j-p/2 : 2,l,l,l,l).(m + a;/2 : 1.0,1.1,1) 
{l + k + m + u;/2:l,0, 0, 1, 1), (1 - j - p/2 : 1, 0, 1, 0, 1) 
(fc + w/2: 1,1,0,1,1) 



T7 3:0;0;0;0;0 
3:0;0;0;0;0 



(l-i-p/2: 1,1,0,1,0) 
where x = p'^/pl. 



/y> /y> /y> nf 



(37) 



4. CONCLUSION 

The technique of Feynman parametrization can of course be used to perform 
loop integrals in different gauges but it is very difficult to perform the parametric 
integrals for arbitrary exponents of propagators. Not so with NDIM, carry loop 
integrals out with particular exponents is as easy as dealing with arbitrary ones - 
besides, one can come across with singularities which depend on them and not on 
dimension D - this fact is very important when we are studying light-cone gauge 
Feynman integrals, because one could have to handle products like ((j^)" [{q — p)^] , 
being a and b negative. NDIM can calculate all of them simultaneously, but if one 
chooses partial fractioning tricks then he/she will be forced to carry out each inte- 
gral separately. Besides usual covariant integrals and the trickier light-cone gauge 
ones, NDIM was probed in the Coulomb gauge, where a procedure introduced 
by Leibbrandt and co-workers - called split dimensional regularization is needed in 
order to render the energy integrals well-defined. 

In this paper, we studied Feynman loop integrals pertaining to three outstanding 
gauges: the usual, and more popular, covariant Feynman gauge and two of the 
trickiest non-covariant gauges, the light-cone and the Coulomb ones. Our results 
are given in terms of hypergeometric functions and in the dimensional regularization 
context. 
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